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Abstract. Using a high-statistics sample of neutral current elastic neutrino interactions, MiniBooNE measured the flux- 
averaged neutral current elastic differential cross-section on mineral oil (CH2). Using the latter, a test of MC with different 
values of the axial vector mass has been performed. Also, a possibility of using a sample of neutral current elastic proton- 
enriched events above Cherenkov threshold to measure the ratio vp vp/vN vN is discussed. This ratio is sensitive to 
the strange quark contribution to the nucleon spin, A^. 
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NEUTRAL CURRENT ELASTIC INTERACTIONS IN MINIBOONE 

Neutrino-nucleus neutral current elastic (NCE) scattering, vA^ vN (N = p, n), is a fundamental probe of the nucleon. 
This interaction, in principle, is sensitive to the isoscalar part of the weak neutral current axial form factor, especially 
when it is possible to separate protons from neutrons. In MiniBooNE, NCE events account for about 18% of all 
neutrino interactions. It seems virtually impossible to measure these interactions in a Cherenkov detector, because 
most of the protons have an energy below Cherenkov threshold. Also, there is not enough light for tagging the neutron 
capture to separate neutrons from protons. However, the use of mineral oil as a medium in the MiniBooNE detector jll] 
allows for a measurement of these interactions even in a region below Cherenkov threshold, because of the presence 
of fluors in the mineral oil. Hence, charged particles going through the medium produce scintillation and fluorescence 
light in addition to any Cherenkov light. 

There was a preliminary measurement of the MiniBooNE NCE differential cross-section |2], where a reduced 
MiniBooNE data set with 6.57 x lO'^ protons on target (POT) was used, at the time restricted by the oscillation 
blindness requirement. Since then, MiniBooNE has improved the analysis, as well as used the full data set with 
6.46 X lO^o POT for the latest NCE results. The NCE flux -averaged differential cross-section has been measured and 
is reported. It has also been used for the tests of MC with different values of the axial vector mass. Ma- 

NCE scattering may occur off either protons or neutrons. In the detector, such an interaction is observed via 
a recoil nucleon. As opposed to protons, neutrons are neutral particles and do not produce light by themselves. 
They are normally observed in the detector via a charge exchange process np np that follows the initial neutrino 
interaction. For this reason, in MiniBooNE it is impossible to distinguish between protons and neutrons for protons 
below Cherenkov threshold in mineral oil (i.e., for proton kinetic energy T ^ 350 MeV). 

Because of the impossibility of measuring the outgoing neutrino in such interactions, a stationary target nucleon is 
assumed in the momentum transfered squared (Q^) reconstruction. By measuring the outgoing nucleon kinematics, 
the reconstructed may be calculated as: 

Q^^lmNT, (1) 

where niN is the nucleon mass and T is the reconstructed nucleon kinetic energy. 

The signal consists of 3 different processes: vA^ — » vA^ scattering off bound protons and neutrons inside carbon 
atoms and off hydrogen atoms. Thus, MiniBooNE measures the NCE flux-averaged differential cross-section, which 
also averages over the signal processes. The weight for each signal process in the total NCE cross-section is given in 
EqM- 

Neutrino interactions within the detector are simulated with the v3 Nuance event generator iH, where the relativistic 
Fermi gas model of Smith and Moniz 101 is used to describe NCE interactions. For resonant pion production the Rein 
and Sehgal model iSfl is used. In the few GeV range, such interactions are dominated by the A(1232) resonance. 



TABLE 1. Cross-section parameters and their 1-C7 uncertainties used in the MiniBooNE MC. 



Parameter 



Value 



Uncertainty Description 



dels 

M_A Ipi 
M_A Npi 
coh 
deh-ad 
dis 

EloSF 
M_A Coh 
Res PiO 
M_A QE H 



M_AQE 

E_B 

P-F 



1.2341 GeV 
34.0 MeV 
220.0 MeV 
0.0 

1.1 GeV 
1.3 GeV 
1.302 
1.00 
1.00 
1.0220 
1.030 GeV 
1.00 

1.13 GeV 



0.077 GeV 
9.0 MeV 
30.0 MeV 
0.1 

0.27 GeV 
0.52 GeV 
0.14 
0.12 
0.25 
0.0205 
0.275 GeV 
0.14 

0.10 GeV 



M/i for QE events on carbon 
Binding energy for carbon 
Fermi momentum for carbon 
As, the axial vector isoscalar term 

for CC and NC Ipi events 
Ma for CC and NC multipion events 
Scale factor for NC coherent piO events 
Scale factor for NC and CC delta radiative events 
Scale factor for deep inelastic scattering events 
Pauli blocking K 

Ma for CClpi coherent events (not NCpiO coherent) 
Scale factor for NC resonant piO events 
Ma for QE events on Hydrogen 



although, contributions from higher mass resonances are also included. The values of the parameters assumed in the 
neutrino cross-section model and their uncertainty are shown in Table [T] 

Intranuclear final state interactions (FSI) inside '^C are also interpreted in Nuance, where the scattered hadrons are 
propagated through the nucleus, which is simulated based on models of nuclear density and Fermi momentum. At each 
step, the MC determines whether the hadron interacts with the target nucleus and if so, generates the interaction and the 
resulting new final state particles. Particle propagation in the detector is a GEANT3-based g] MC with GCALOR 01 
hadronic interactions simulating the detector response to particles produced in the neutrino interactions. Due to the 
FSI, a NCE event may have more than 1 final state particle. Also, a NC pion event may have no pion in the final state 
as the pion can get absorbed within the carbon nucleus. According to Nuance, for NCE scattering on a target nucleon 
inside '^C, the probability of having only 1 outgoing nucleon is ^ 74%. Also, the probability that the outgoing pion 
experiences either absorption or charge exchange is ^ 20%. 

Neutrino events that happen outside of the detector in the dirt or in the detector material (referred to as "dirt" 
background henceforth) are simulated the same way as in-tank events but with a cross-section reweighed according to 
the density of the material relative to the density of the mineral oil. 



EVENT RECONSTRUCTION AND SELECTION. DIRT BACKGROUND CONSTRAINT. 



One of the main analysis improvements is the NCE event reconstruction, which assumes that an event is a single 
proton. The fitter looks for a set of parameters (position, time, direction and energy) which maximizes the event 
likelihood using the charge and time information from all photomultiplier tubes (PMTs). The event is assumed to be 
point-like with the Cherenkov and scintillation profiles for protons measured from the MC. The position resolution is 
^ 0.75 m, the energy resolution ^ 25%, and for protons above Cherenkov threshold the angular resolution is ~ 10". 
The reconstructed nucleon kinetic energy, T, is proportional to the total charge on all PMTs, and thus it is a measure 
of the sum of the kinetic energies of all final state nucleons that are produced in the interaction. 

The following set of cuts are applied to the full MiniBooNE data set to select the NCE sample: (1) only 1 subevent 
to ensure no decaying particles {ji or K decay), (2) veto PMTs < 6 to remove events exiting or entering the detector, (3) 
tank PMTs > 25 to use only reconstructable events, (4) beam time window cut, (5) a cut on the reconstructed proton 
energy T < 650 MeV (above which the NCE signal to background ratio decreases significantly), (6) a cut on the time 
log-likelihood ratio between proton and electron event hypotheses ln(^e/^p) < 0.42 to eliminate beam-unrelated 
electrons from cosmic ray muon decays (From FiglTJa) one can see that electrons have more prompt light fraction 
than protons), (7) a fiducial volume cut, R < 4.2 m if T < 200 MeV and R < 5 m otherwise. The reason to have a 
tighter fiducial volume at low energies is to reduce dirt backgrounds. 

A total of 94,531 events pass the NCE cuts in the sample coiTesponding to 6.46 x 10^*^ POT. The efficiency of the 
cuts is estimated to be 26% (a large part of the reduction is due to the fiducial volume cut). The fraction of NCE events 
in the sample is 65%. The remaining 35% of events in the sample are backgrounds: 15% irreducible backgrounds, 
10% dirt events, and 10% other backgrounds (of which only 0.5% are beam unrelated). The reconstructed nucleon 
kinetic energy spectrum for the selected NCE event sample is shown in FiglUb). 




FIGURE 1. (a): Prompt light fraction (hits with —5 ns < tcorr < 5 ns) for the beam unrelated data and the NCE MC under electron 
hypothesis reconstruction, (b): Reconstructed nucleon kinetic energy spectra for the data and different MC channels after the NCE 
selection cuts applied, but using a uniform fiducial volume cut of R < 4.2 m. Distributions are absolutely (POT) normalized. 



Irreducible backgrounds are those neutral current n production channels that have no TT in the final state, as the 
n is absorbed in the nucleus through a FSI process. In this case, the final state particles for these events are just 
nucleons, identical to the final state particles produced in NCE events and hence are "irreducible". These backgrounds 
are important at intermediate energies, 200 MeV< T < 500 MeV. 

The dirt background is an important contribution to the NCE data sample at low energies, mostly below 200 MeV. 
Normally, this background is due to nucleons (mainly neutrons), which are produced in neutrino interactions outside 
of the detector, penetrating into the detector without firing enough veto PMTs. Dirt events are difficult to simulate in 
the MC because they happen outside of the detector in the soil, whose properties have not been studied in detail. It is 
also hard to correctly describe the particle propagation through a set of different materials (through the soil, detector 
parts and into the detector media). However, it is possible to directly measure these backgrounds in the data, as they 
have different kinematic observables than in-tank events, such as the reconstructed radius R (dirt events are mostly 
in outer regions), the reconstructed Z (dirt events are mostly upstream with Z < 0) and the reconstructed energy (dirt 
events are at T < 200 MeV). The dirt energy spectrum is measured by fitting the MC in-tank and dirt templates to 
the data in Z, R, and energy distributions for the dirt-enriched event samples. All of these measurements agree with 
each other within 10%. The measured number of dirt events is determined to be ^ 30% lower than the original MC 
prediction. 



FINAL NEUTRAL CURRENT ELASTIC DIFFERENTIAL CROSS-SECTION RESULT. 

The obtained reconstructed energy spectrum is then corrected for limited detector resolution and detector efficiency 
effects that distort the original spectrum. This is the so-called unfolding problem i^, for which a method based on 
Bayes' theorem is used. This method gives a well-behaved but biased solution which depends on the original MC 
energy spectrum. The eiTor due to the bias is estimated by iterating the unfolding procedure, where the new MC energy 
spectrum is replaced by the unfolded energy spectrum. 

Possible systematic uncertainties and their contribution to the total error have been studied and are shown in Table|2] 
With this, the NCE flux-averaged differential cross-section is shown in Fig|2] as a function of Q^. The predicted 
distribution of the irreducible backgrounds, which have been subtracted along with the rest of the backgrounds, is also 
shown in the figure. 

The average chemical formula for mineral oil is CH2. Thus, the vA^ vA^ cross-section in terms of the NCE 
cross-section on different types of bound and free nucleons is expressed as: 

-J^ - ^C,,^vp,h{Q ) +jC,p^vp.dQ ) ^Q2 +7Cv„^v«r(e ) ^Q2 ' (2) 

where the first term is the NCE cross-section on free protons taken per free proton, etc. The efficiency correction 
functions Cf are due to different efficiencies for the signal processes. They are defined as the ratio of the efficiency 



TABLE 2. The individual error contributions to the total error for the reconstructed 
energy spectrum and true energy spectrum (or the NCEL cross-section error). 



Error 


Value, % 


Statistical 


3.1 


Discriminator threshold 


1.1 


QT PMT response 


2.0 


Dirt 


1.4 


POT 


1.7 


7t" yield 


0.1 


Hadronic 


0.2 


Cross-section 


3.4 


Horn uncertainties 


4.7 


A^o production 


0.1 


production 


0.5 


7t^ production 


0.3 


7t^ production 


4.1 


Optical model 


12.3 


Unfolding 


7.5 


Total 


18.9 
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FIGURE 2. MiniBooNE vA' vN flux-averaged differential cross-section on CH2 along with the predicted irreducible back- 
ground cross-section, that has been subtracted out. 



for a particular type of signal events to the average efficiency for all signal events. In principle, one needs to know the 
efficiency corrections Cyp H, Cyp^c ™d Cy„,c, which are shown in Fig|3]!a). But for some quick calculations one can 
roughly assume Cyp,// = Cvp,c = Cv«,C = 1- 

The NCE cross-section is quite sensitive to the axial form factor In fact, at low Q^, da/dQ^ ^ G\{Q^). The axial 
form factor is usually taken in the dipole approximation: 

GA(e') = Q.5{TgA+^s)/{\ + Q^/Mlf (3) 

where gA — ^ -267 1 , Ma is the axial vector mass, the + sign is for protons, and the — sign is for neutrons. 

Inserting Eq.© for protons and neutrons into Eq.(|2|i, one can see that the MiniBooNE vA^ vN NCE cross-section 
is not sensitive to As, as the linear term in As nearly cancels, while the quadratic term in As remains, but it is small if 
I As I <C gA- However, one can still use these data to probe Ma- 

There seems to be some recent interest in the value of Ma- Before MiniBooNE, Ma has been measured by 
mostly bubble chamber deuterium-based neutrino experiments, which found Ma — 1 .02 GeV [|10<1. However, recent 
measurements on carbon, oxygen and iron targets by MiniBooNE [|l 1|] . K2K lfT2ilT3ll and MINOS |[3 have determined 




FIGURE 3. (a): Efficiency corrections, Cf, as a function of the nucleon kinetic energy for different NCE processes, (b): -tests 
on the NCE reconstructed energy for values of 1.35, 1.23, and 1.02 GeV. The values are 27.1, 29.2 and 41.3 for 49 DOF, 
respectively. The distributions are absolutely normalized. 



it to have ^ 20 — 30% larger values, using the CCQE channel. In the MiniBooNE MC it is possible to generate the 
reconstructed energy spectra for NCE samples that correspond to any value of Ma- Here we test how different values of 
Ma agree with the data assuming absolutely normalized distributions of nucleon energy. To calculate the between 
data and MC, the full error matrix has been used which takes into account the variations of MC energy spectrum due 
to various MC uncertainties (such as flux, cross-section, beam etc). The reconstructed energy spectrum for the data 
and MC for M^ and K values that we are interested in are shown in FigOb). The higher values of M^ seem to describe 
the MiniBooNE NCE data better 



HIGH ENERGY PROTON-ENRICHED SAMPLE. A5 SENSITIVITY. 

As mentioned before, the Ai sensitivity of the NCE sample comes down to the possibility of distinguishing proton 
from neutron events. In particular, a ratio of vp vp to vn vn is most sensitive to As in addition to reducing the 
systematic error 

In MiniBooNE, it is virtually impossible to distinguish protons from neutrons, because neutrons are usually only 
registered after a strong interaction with a proton. However, it is possible to select a special class of NCE protons, 
namely those with only 1 proton in the final state and energies above the proton Cherenkov threshold. These events 
can be used for a As measurement relative to the total measured vA^ vA^ cross-section. 

It should be noted that events with multiple nucleons in the final state (NCE neutrons and irreducible backgrounds) 
have multiple protons produced in the reaction. At kinetic energies above 350 MeV, single proton events produce on 
average a higher Cherenkov light fraction than multiple proton events. In addition, these classes differ in the kinematics 
of the outgoing nucleon, with single proton events being more forward-going. 

In order to reduce systematic uncertainties, a ratio vp ^ Vp/vN vN as a function of reconstructed nucleon 
kinetic energy from 350 MeV to 800 MeV is used. The denominator of the ratio are events with NCE selection cuts 
used in the cross-section measurement, but the energy cut (5) has been replaced with 350 MeV< T < 800 MeV, and 
an additional cut is implemented, ln(^^/^;;,) (a log likelihood ratio between muon and proton hypotheses) in order to 
reduce muon-like backgrounds that dominate at high visible energy region. With these cuts, the vA^ vN data sample 
has 24,004 events with predicted channel fractions: 45% NCE, 26% irreducible backgrounds, 3% dirt backgrounds 
and 25% other backgrounds. 

For the numerator, two more cuts are applied in addition to the ones used for the denominator. The first cut is based 
on the prompt light fraction. It reduces multiple proton events, which have less Cherenkov light than single proton 
events. After this cut, looking at the Op distribution (the angle between the reconstructed nucleon direction and the 
beam direction), shown in FigHJa), one can see that most of the single proton events are forward going. Thus, the last 
cut applied for the NCE p sample is Op < 60". The final vp Vp data sample has 7,616 events with predicted channel 
fractions: 55% NCEp, 10% NCE n, 14% irreducible backgrounds, 1% dirt backgrounds, and 19% other backgrounds. 




FIGURE 4. (a): Op, the angle between the nucleon reconstructed direction and the beam direction for data and MC after NCE 
and p/ fi cuts are applied. All distributions are absolutely normalized, (b): The ratio vp vp/vN —> vN as a function of the 
reconstructed energy for MC with Ai- = with total error. Also MC with As = —0.5 and 0.5 are shown. 



The ratio of vp ^ Vp to vA' vN for the MC with A.s- = and the total uncertainties along with the MC with values 
of As — —0.5 and 0.5 is shown in FigHJb). 

This will be the first attempt at a Ai measurement using this ratio. The systematic errors are quite large, mostly 
due to large uncertainties in the optical model of the mineral oil. MiniBooNE has sensitivity to Ai for proton energies 
above Cherenkov threshold for protons, where the contribution from irreducible backgrounds is significant. In order to 
have a good sensitivity for As, future experiments need to have good proton/neutron identification (possibly through 
neutron capture tagging) and extend down T < 200 MeV region, where the contribution from irreducible backgrounds 
is less. 

In summary, MiniBooNE has used a high-statistics sample of NCE interactions to measure the total NCE cross- 
section. Using absolutely normalized distributions of the reconstructed energy for the NCE sample, tests for several 
Ma and k values have been performed. The MC with higher values of Ma, 1.23 GeV and 1.35 GeV gives a better 
X^, than 1.02 GeV, which is in agreement with the shape normalized fits of the CCQE channel obtained by recent 
experiments. For energies above Cherenkov threshold, a sample of NCE proton-enriched protons was obtained. These 
events allow a measurement of the ratio of vp ^ vp to vA^ vN and is sensitive to As. The approximate sensitivity 
plot for the MiniBooNE Ai extraction was shown. 
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